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Abstract 

The increasing number of objects orbiting the Earth justifies the great attention and interest in the ob- 
servation, spacecraft protection and colhsion avoidance. These studies involve different disturbances and 
(~| ■ resonances in the orbital motions of these objects distributed by the distinct altitudes. In this work, the 

TLE (Two-Line Elements) of the NORAD are studied observing the resonant period of the objects orbiting 
^ ■ the Earth and the main resonance in the LEO region. The time behavior of the semi-major axis, eccentricity 

O I and inclination of some space debris are studied. Possible irregular motions are observed by the frequency 

■ analysis and by the presence of different resonant angles describing the orbital dynamics of these objects. 
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1. Introduction 
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pi ^ The objects orbiting the Earth are classified, basically, in Low Earth Orbit (LEO), Medium Earth Orbit 

(MEO) and Geostationary Orbit (GEO). Most of the objects are found in the LEO region because this region 
has a big quantity of space debris. Considering approximately 10000 cataloged objects around the Earth, 
one can verify the distribution of the objects as: 27% of operational spacecraft, 22% of old spacecraft, 41% of 
miscellaneous fragments, 17% of rocket bodies and about 13% of mission-related objects. The uncatalogued 
objects larger than 1 cm are estimated in some value between 50000 and 600000 [l[[2|- 

Currently, the orbital motions of the cataloged objects can be analyzed using the 2-line element set of 
the NORAD (North American Defense) fsl] . The TLE are composed by seven parameters and epoch. These 
data can be compared, for example, with the model of the orbit propagator situated in the artificial satellite. 
A similar study is done for the Brazilian satellite CBERS-1 in cooperation with China. In this case, orbital 
perturbations due to geopotential, atmospheric drag, solar radiation pressure, gravitational effects of the 
Sun and the Moon are considered in the numerical integration of the orbit and the results are compared 
with the TLE data 

In the last years, the LEO region have been studied about the space debris mitigation due to the increasing 
number of this kind of object through the years. These aspects englobe the observation, spacecraft protection 
and collision avoidance 0, @|- The space debris are composed of aluminum from spacecraft structures, 
alumina from solid rocket motor exhausts, zinc and titanium oxides from thermal control coatings and their 
size ranges from several meters to a fraction of a micrometer in diameter [9]. 

The space between the Earth and the Moon has several artificial satellites and distinct objects in some 
resonance. Synchronous satellites in circular or elliptical orbits have been extensively studied in literature, 
due to the study of resonant orbits characterizing the dynamics of these satellites since the 60's Lia li, il 

n [11, [n [n, M, M, Ei, s [m, m m, h m, m, m,^^- 
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In this work, objects in resonant orbital motions are investigated. This search involves all cataloged 
objects in the TLE files of the NORAD. Figures show regions with objects in deep resonance. The time 
behavior of the semi-major axis, eccentricity, inclination and the frequency analysis corresponding to the 
orbital motions of some space debris are studied. Possible irregular motions are observed by the presence of 
different resonant angles describing the orbital dynamics of these objects. 



2. Resonant objects orbiting the Earth 

In this section, the TLE data are used to verify objects in resonant orbital motions and what resonance 
has the majority of objects [3]. 

The present distribution of objects by the value of the mean motion n indicates the commensurability 
between the frequencies of the mean motion of the object and the Earth's rotation motion. See the histogram 
of the mean motion in Fig. [TJ 
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Figure 1: Histogram of the mean motion of the cataloged objects. 

In Fig. [T]it is verified that most of the objects are in the region 13 < n{rev / day) < 15. 
To study the resonant objects using the TLE data, a criterium is established in the resonant period 
Pres, by the condition Pres > 300 days. Pres is obtained by the relation. 



Pres 



and 4>impq is calculated from (l^ . 



27r 

4^lmpq 



(1) 



<j)i^pq{M,uj, n, e)^{l-2p + q)M + 2p)uj + m{n - e - A,„) + (/ - m)- . (2) 

with a, e, /, fi, uj, M are the classical keplerian elements: a is the semi-major axis, e is the eccentricity, 
I is the inclination of the orbit plane with the equator, ft is the longitude of the ascending node, w is the 
argument of pericentre and M is the mean anomaly, respectively; & is the Greenwich sidereal time and Xim 
is the corresponding reference longitude along the equator. So, (pimpg is defined as 
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'Pi^pg = {l-2p + q)M 2p)uj + m{h - 6) . (3) 

Substituting fc = Z — 2p in ([3]) , 



(t^hmq = {k + q)M + kuj + m{ri — 9) 
The terms w, tl and M can be written as Id, Isdl. 



(4) 
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Ce is the Earth mean equatorial radius, ae=6378.140 km, J2 is the second zonal harmonic, J2 = 1, 0826 x 10 ^. 
The term in rad/day is 



e w 1.00 2 7 3 7 9 9 26 x 27r. 



(6) 



In order to use orbital elements compatible with the way in which Two-Line Elements were generated, 
some corrections are done in the mean motion of the TLE data. Considering as rii the mean motion of the 
2-line, the semi-major axis ai is calculated ^]. 
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where n is the Earth gravitational parameter, /i=3. 986009 x 10^^ w? j . Using ai, the parameter 8\ is 
calculated by the Eq. dH) 



_ 3 (3cog^(/)-l) 
'"4 'a? (i_e2)3/2 ' 



(8) 



Now, the new semi-major axis Oo used in the calculations of the resonant period is defined using b\ from 
the Eq. (El) [i. 
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and the new mean motion Uo used in the calculations is found considering the semi-major axis corrected ao 
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Now, using the corrections previously shown, a TLE file is studied with the purpose to investigate what 
resonance (defined by the commensurability between the mean motion of the object and the Earth's rotation 
angular velocity) has most of the objects orbiting the Earth. After this analysis some space debris in deep 
resonance are investigated in the third section. 

The file analyzed is "a^Wata_2011_045" from the website of the Space Track [3], and it corresponds to 
february 2011. In Tab. 1, the number of data and the number of different objects are specified. 



Table 1: 2-line data of objects orbiting the Earth. 

File Number of data Number of different objects 

alldata 2011 045 15360 9745 



The simulation identified objects with resonant period greater than 300 days. Several values of the 
coefficients, fc, q and m are considered in the Eq. ([2]) producing different resonant angles to be analyzed by 
the Eq. dl]). See Tab. 2 showing details about the results of the considered simulation. 



Table 2: Results of the simulation: objects in resonant orbital motions 



Number of data 


Number of different objects 


coefficient k 


coefficient q coefficient m 


3180 


1276 


-50 < A: < 50 


-5<q<5 l<m<50 



Note that Tab. 2 shows only the data satisfying the established criterium Pres > 300 days. The results 
show the resonant angles and resonant periods which compose the orbital motions of the resonant objects. 

Figs. [5] and [3] show the semi-major axis versus eccentricity and the semi-major axis versus inclination, 
respectively, using only the data of the related objects in Tab. 2. 



objects 



-% leeee 



: Hf 



IBBBB '^ 



+ +M#HjmH,J- -H- _+h+ 



e e.i B.2 B.3 B.4 a. 5 e.6 b.7 B.a e.9 

Eccentridt)' 



Figure 2: Semi-major axis versus eccentricity of resonant objects satisfying the condition of the Pres > 300 days. 
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Figure 3: Semi-major axis versus inclination of resonant objects satisfying tlie condition of the Pres > 300 days. 



The objects with Pres > 300 days can be visualized around the Earth in three different regions by the 
value of the semi-major axis, a < ISOOOfcm, 15000A:to < a < 40000A:to and a > 40000/cm, see Fig. U) 




Figure 4: Objects around the Earth satisfying the condition of the resonant period greater than 300 days. 
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Tab. [3] shows values of the coefficient m and the number of objects. It is possible to verify that most of the 
orbital motions of objects are related with the coefficient m = 14 considering the value of the semi-major 
axis up to ISOOOfcm. 

Table 3; Number of objects by the number of the coefficient m considering the value of the semi-major axis up to iSOOOfcm 



TO 


Number of objects 


TO Number of objects 


7 


2 


30 


31 


8 


1 


31 


13 


9 


3 


32 


5 


10 


1 


33 


2 


11 


3 


34 


1 


12 


20 


35 


4 


13 


31 


36 


5 


14 


351 


37 


42 


15 


65 


38 


14 


16 


4 


39 


14 


18 


3 


40 


24 


20 


1 


41 


186 


21 


2 


42 


116 


22 


2 


43 


114 


23 


2 


44 


52 


24 


10 


45 


25 


25 


93 


46 


8 


26 


19 


47 


4 


27 


66 


48 


6 


28 


179 


49 


22 


29 


97 


50 


33 



Comparing the number of objects satisfying the condition of Pres > 300 days in the different regions, 
a < 15000fcTO and a > 15000fcTO, it is verified that about 62.37% of the resonant objects has a < 15000A:to. 
See this information in Tabs. |4]and[5l 



Table 4: Number of objects satisfying the condition Pres > 300 days 


a < 15000fcTO 


15000fcTO <a< 40000fcTO a > 40000fcTO 


1276 


331 439 


Table 5: Percentaj 


;e of objects satisfying the condition Pres > 300 days 


a < 15000fcTO 


15000fcTO < a < 40000fcTO a > 40000fcTO 


62.37 % 


16.18 % 21.46 % 



These studies allow to investigate the real influence of the resonance effect in the orbital dynamics of the 
artificial satellites and space debris. The number of resonant objects in comparison with the total number 
of objects in the TLE data shows the great influence of the commensurability between the mean motion of 
the object and the Earth's rotation angular velocity in its orbits. In this way, a more detailed study about 
the resonant period and the resonant angles is necessary. 

In the next section, the orbital motions of some space debris in deep resonance are studied. 
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3. Study of objects in deep resonance 



Considering the results of the simulations shown in the second section, some cataloged space debris are 
studied with respect to the time behavior of the semi-major axis, eccentricity, inclination, resonant period, 
resonant angle and the frequency analysis of the orbital elements. These data are analyzed in this section 
observing the possible regular or irregular orbital motions. 

Fig. O shows the semi-major axis versus resonant period of objects satisfying the condition Pres > 300 
days. The objects are distributed by the value of semi-major axis in three different regions: 1) a < 15000fcm, 
2) ISOOOfcm <a< 40000fcrn and 3) a > 40000fcm. 
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Figure 5: Semi-major axis versus resonant period of objects satisfying the condition Pres > 300 days. 



Analyzing Fig. [51 it is possible to verify some objects in deep resonance, in other words, the orbital 
motions of these objects are influenced by the resonance for a long time. In order to study the orbital 
motions of these objects, four cataloged space debris are analyzed. From the TLE data, objects are identified 
by the numbers 325, 546, 2986 and 4855. These objects were chosen because they show resonant period 
Pres > 10000 days. 

Figs. to [T71 show the time behavior of the semi-major axis, eccentricity and inclination of the objects 
325, 546, 2986 and 4855, and the frequency analysis of these orbital elements. 

The frequency analysis is used as a tool to verify the chaotic orbits. Several authors. 



34, 35, 3' 



[311,131132 

also use this method to study the possible regular or irregular orbits in different dynamical systems. A brief 
description of the frequency analysis involving the FFT (Fast Fourier Transforms), described with more 



details in |31l. l32l. l36fl. is presented in what follows. 

Generally, the frequency analysis is applied in the output of the numerical integration in the study of 
dynamical systems. However, in the present work, this analysis is used in the real data of the orbital motions 
of the objects 325, 546, 2986 and 4855, from the TLE data of the NORAD [3]. 



For regular motions, the orbital elements oe{t) show a dependence on time as follows 3ll. l32l. l3' 



oe(t) =J2^^' 



27rihft 



(11) 
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where Ah represent the amphtudes, h G and f is a frequency vector. The components of f compose 
the fundamental frequencies of motion, and the spectral decomposition of the orbital motion is obtained 
from the Fourier transform when the independent frequencies are constant in the course of time [36|. It is 
possible to observe that f depends on the kind of trajectory (3l| . 

The same numerical procedure is done for irregular and regular motions. The difference in behavior of 
these orbital motions is used to identify the two kinds of trajectories. The power spectrum, produced by 
the Fourier transform, of regular motion can be distinguished because it generally has a small number of 
frequency components. The irregular trajectories are not conditionally periodic and they do not compose 
an invariant tori. In this way, the Fourier transform of an orbital element, for example, is not a sum over 
Dirac 5- functions, and consequently the power spectrum is not discrete for irregular motions [31., ,36J. 
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Figure 6: Orbital motion of object 325 corresponding to January - September, 2011; a) Time behavior of the semi-major axis 
and b) Power spectrum of the semi-major axis. 
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Figure 7: Orbital motion of object 325 corresponding to January - September, 2011; a) Time behavior of the eccentricity and 
b) Power spectrum of the eccentricity. 
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Figure 8: Orbital motion of object 325 corresponding to January - September, 2011; a) Time behavior of the inclination and 
b) Power spectrum of the inclination. 




100 150 
Time (days), Jan -Sep, 2011 




Figure 9: Orbital motion of object 546 corresponding to January - September, 2011; a) Time behavior of the semi-major axis 
and b) Power spectrum of the semi-major axis. 
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Figure 10: Orbital motion of object 546 corresponding to January - September, 2011; a) Time behavior of the eccentricity and 
b) Power spectrum of the eccentricity. 
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Figure 11: Orbital motion of object 546 corresponding to January - September, 2011; a) Time behavior of the inclination and 
b) Power spectrum of the inclination. 
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Figure 12: Orbital motion of object 2986 corresponding to January - September, 2011; a) Time behavior of the semi-major 
axis and b) Power spectrum of the semi-major axis. 
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Figure 13: Orbital motion of object 2986 corresponding to January - September, 2011; a) Time behavior of the eccentricity 
and b) Power spectrum of the eccentricity. 



11 





100 150 
Time (days), Jan -Sep, 2011 



0.3 0.4 
Frequency 



Figure 14: Orbital motion of object 2986 corresponding to January - September, 2011; a) Time behavior of the inclination and 
b) Power spectrum of the inclination. 
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Figure 15: Orbital motion of object 4855 corresponding to January - September, 2011; a) Time behavior of the semi-major 
axis and b) Power spectrum of the semi-major axis. 
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Figure 16: Orbital motion of object 4855 corresponding to January - September, 2011; a) Time behavior of the eccentricity 
and b) Power spectrum of the eccentricity. 
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Figure 17: Orbital motion of object 4855 corresponding to January - September, 2011; a) Time behavior of the inclination and 
b) Power spectrum of the inclination. 



Observing the time behavior of the orbital elements of the objects 325, 546, 2986 and 4855 in Figs. [6] 
to [T71 one can verify possible regular and irregular motions in the trajectories of these space debris. The 
time behavior of the inclination of the orbital motions of the objects 325, 546 and 2986 show irregularities. 
Analyzing object 325, note that in 200 days, a fast increase in the inclination occurs but, this variation is 
about 0.01° and it may be related with some disturbance added to the motion. 

The power spectrum of the orbital elements are not discrete and they have a big number of frequency 
components in Figs. [6] to [171 Usually, in the studies of the dynamical systems with the purpose to identify 
chaos, the frequency analysis is considered for a long time. But, in Low Earth Orbits, a big number of 
cataloged objects have their motions influenced by different resonant periods and unknown objects can 
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collide with each other providing more unpredictability in their motions. In this case, the frequency analysis 
can be more appropriated for a short time. 

So its also important to observe if the orbital motion is influenced by different resonant angles. In other 
words, it is verified if the orbital motion is due to overlap resonances. 

Tab. [5] shows the resonant angles related to the orbital motions of the objects 325, 546, 2986 and 4855 
corresponding to the period January - September, 2011. 

Table 6: Resonant angles ipkmq related to the orbital motions of the objects 325, 546, 2986 and 4855. 



Object 


Coefficient k 


Coefficient m 


Coefficient q 


325 


1 


14 





325 


6 


42 


-3 


325 


4 


28 


-2 


325 


2 


14 


-1 


325 


7 


42 


-4 


325 


5 


28 


-3 


325 


3 


14 


-2 


325 


8 


42 


-5 


325 


6 


28 


-4 


325 


4 


14 


-3 


325 


7 


28 


-5 


546 


3 


43 





546 


4 


43 


-1 


546 


5 


43 


-2 


546 


6 


43 


-3 


546 


7 


43 


-4 


546 


8 


43 


-5 


2986 


-1 


28 


3 


2986 


-1 


42 


4 


2986 





14 


1 


4855 


7 


41 


-4 


4855 


8 


41 


-5 



Observing Tab. [51 one can verify that the orbital dynamics of the objects 2986 and 4855 are influenced 
by some resonant angles, while for the objects 325 and 546, several resonant angles influence their orbits 
simultaneously. 

If the commensurability between the orbital motions of the object and the Planet is deflned by the 
parameter a and by the condition a = {k + q)/m, one can say that the exact 14:1 resonance is deflned by 
the condition a = 1/14. This way, analyzing Tab. [51 it is verified that the motions of the objects 325 and 
2986 are influenced by the exact 14:1 resonance, while the objects 546 and 4855 are influenced by resonant 
angles in the neighborhood of the exact resonance. 

Figures [15] to [H] show the time behavior of the resonant period corresponding to the resonant angles 
presented in Tab. [51 
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Figure 18: Time behavior of the resonant period corresponding to the orbital motion of object 325. 
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Figure 19: Time behavior of the resonant period corresponding to the orbital motion of object 546. 
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Figure 20: Time behavior of the resonant period corresponding to the orbital motion of object 2986. 
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Figure 21; Time behavior of the resonant period corresponding to the orbital motion of object 4855. 
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Analyzing the time behavior of the resonant period in Figs. [18] to [2T] it is verified that the resonant 
angles remain in deep resonance for a few days. The orbital motion of object 546 has resonant angles which 
satisfies the established criterium Pres > 300 days for 70 days and consequently the orbital dynamics of this 
object may be losing the infiuence of the 14:1 resonance. Figures [50] and [^T] corresponding to the objects 
2986 and 4855 show the resonant period for values greater than 10000 days in the first fifty days, after that 
decreasing for values smaller than 500 days in the resonant period and it shows a tendency to not remain 
in deep resonance. The orbital motion of object 325 is more influenced by the 14:1 resonance and this fact 
can be observed by the time behavior of the (j>kmq; Fig- 1221 

To continue the analysis about the irregular orbital motions, the time behavior of the (jjkmq is studied 
verifying if different resonant angles describe the orbital dynamics of these objects at the same moment. 
This fact can indicate irregular orbits requiring a full system with all resonant angles as described in Tab. 
[6] Figures [22] to [25] show the time behavior of the (pkmq corresponding to the resonant angles presented in 
Tab. E 
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Figure 22: Time behavior of <t>km,q corresponding to tlie orbital motion of object 325. 
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Figure 23: Time behavior of <j)kmq corresponding to the orbital motion of object 546. 
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Figure 24: Time behavior of 4'kmn corresponding to the orbital motion of object 2986. 
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Figure 25: Time behavior of ( 



corresponding to the orbital motion of object 4855. 



Analyzing the time behavior of the resonant angles (pkmq in Figs. [22] to [25] it is possible to observe 
irregular orbital motion for object 325. In the beginning, the orbital motion of object 546 is governed by 
distinct resonant angles but, after seventy days, the influence of resonant effects decreases or disappears. 
Objects 325 and 2986 have their orbital motions influenced by the exact 14:1 resonance and they need a full 
system with different resonant angles which compose their motions. Otherwise, the orbital motion of object 
4855 is defined for resonant angles separately, which can be defined as a regular orbit. 

The results and discussions show the complexity in the orbital dynamics of these objects caused by the 
resonance effects. Furthermore, the increasing number of space debris and collisions between them can cause 
big problems for artificial satellites missions. 



4. Conclusions 

In this work, the orbital dynamics of synchronous space debris are studied. From the TLE data of the 
NORAD, objects orbiting the Earth in deep resonance are investigated. 

Analyzing the cataloged objects satisfying the established criterium of the resonant period greater than 
300 days, four space debris (325, 546, 2986 and 4855) are studied observing the irregular characteristics in 
their orbits. The time behavior of the orbital elements and the respective frequency analysis are used to 
verify the orbits of these objects. 

The time behavior of the resonant period of the objects show large amplitudes, and consequently, the 
resonant angles alternate between deep and short resonance. 

Several resonant angles influence, simultaneously, the orbital motions of objects 325 and 2986. The 
time behavior of the inclination and the frequency analysis confirm the irregularities in their orbits. The 
orbital dynamics of object 546 is firstly influenced by several resonant angles in the neighborhood of the 14:1 
resonance but, after some days, the regular characteristic in its orbit seems to dominate. The time behavior 
of the (j)kmq of object 4855 shows resonant angles separately indicating a regular orbit. 

The TLE catalog of the objects orbiting the Earth is an important tool to study the orbital motion of 
the artiflcial satellites and space debris. 
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